VRK1 Regulates Sensitivity to Oxidative Stress by Altering Histone Epigenetic Modifications and the Nuclear Phosphoproteome in Tumor Cells.

The chromatin organization and its dynamic remodeling determine its accessibility and sensitivity to DNA damage oxidative stress, the main source of endogenous DNA damage. We studied the role of the VRK1 chromatin kinase in the response to oxidative stress. which alters the nuclear pattern of histone epigenetic modifications and phosphoproteome pathways. The early effect of oxidative stress on chromatin was studied by determining the levels of 8-oxoG lesions and the alteration of the epigenetic modification of histones. Oxidative stress caused an accumulation of 8-oxoG DNA lesions that were increased by VRK1 depletion, causing a significant accumulation of DNA strand breaks detected by labeling free 3'-DNA ends. In addition, oxidative stress altered the pattern of chromatin epigenetic marks and the nuclear phosphoproteome pathways that were impaired by VRK1 depletion. Oxidative stress induced the acetylation of H4K16ac and H3K9 and the loss of H3K4me3. The depletion of VRK1 altered all these modifications induced by oxidative stress and resulted in losses of H4K16ac and H3K9ac and increases in the H3K9me3 and H3K4me3 levels. All these changes were induced by the oxidative stress in the epigenetic pattern of histones and impaired by VRK1 depletion, indicating that VRK1 plays a major role in the functional reorganization of chromatin in the response to oxidative stress. The analysis of the nuclear phosphoproteome in response to oxidative stress detected an enrichment of the phosphorylated proteins associated with the chromosome organization and chromatin remodeling pathways, which were significantly decreased by VRK1 depletion. VRK1 depletion alters the histone epigenetic pattern and nuclear phosphoproteome pathways in response to oxidative stress. The enzymes performing post-translational epigenetic modifications are potential targets in synthetic lethality strategies for cancer therapies.

Kinase activities in pancreatic ductal adenocarcinoma with prognostic and therapeutic avenues.

Pancreatic ductal adenocarcinoma (PDAC) is a devastating disease with a limited number of known driver mutations but considerable cancer cell heterogeneity. Phosphoproteomics provides a direct read-out of aberrant signaling and the resultant clinically relevant phenotype. Mass spectrometry (MS)-based proteomics and phosphoproteomics were applied to 42 PDAC tumors. Data encompassed over 19 936 phosphoserine or phosphothreonine (pS/T; in 5412 phosphoproteins) and 1208 phosphotyrosine (pY; in 501 phosphoproteins) sites and a total of 3756 proteins. Proteome data identified three distinct subtypes with tumor intrinsic and stromal features. Subsequently, three phospho-subtypes were apparent: two tumor intrinsic (Phos1/2) and one stromal (Phos3), resembling known PDAC molecular subtypes. Kinase activity was analyzed by the Integrative iNferred Kinase Activity (INKA) scoring. Phospho-subtypes displayed differential phosphorylation signals and kinase activity, such as FGR and GSK3 activation in Phos1, SRC kinase family and EPHA2 in Phos2, and EGFR, INSR, MET, ABL1, HIPK1, JAK, and PRKCD in Phos3. Kinase activity analysis of an external PDAC cohort supported our findings and underscored the importance of PI3K/AKT and ERK pathways, among others. Interestingly, unfavorable patient prognosis correlated with higher RTK, PAK2, STK10, and CDK7 activity and high proliferation, whereas long survival was associated with MYLK and PTK6 activity, which was previously unknown. Subtype-associated activity profiles can guide therapeutic combination approaches in tumor and stroma-enriched tissues, and emphasize the critical role of parallel signaling pathways. In addition, kinase activity profiling identifies potential disease markers with prognostic significance.

Protein kinase D drives the secretion of invasion mediators in triple-negative breast cancer cell lines.

The protein kinase D (PKD) family members regulate the fission of cargo vesicles at the Golgi complex and play a pro-oncogenic role in triple-negative breast cancer (TNBC). Whether PKD facilitates the secretion of tumor-promoting factors in TNBC, however, is still unknown. Using the pharmacological inhibition of PKD activity and siRNA-mediated depletion of PKD2 and PKD3, we identified the PKD-dependent secretome of the TNBC cell lines MDA-MB-231 and MDA-MB-468. Mass spectrometry-based proteomics and antibody-based assays revealed a significant downregulation of extracellular matrix related proteins and pro-invasive factors such as LIF, MMP-1, MMP-13, IL-11, M-CSF and GM-CSF in PKD-perturbed cells. Notably, secretion of these proteins in MDA-MB-231 cells was predominantly controlled by PKD2 and enhanced spheroid invasion. Consistently, PKD-dependent secretion of pro-invasive factors was more pronounced in metastatic TNBC cell lines. Our study thus uncovers a novel role of PKD2 in releasing a pro-invasive secretome.

Loss of VRK1 alters the nuclear phosphoproteome in the DNA damage response to doxorubicin.

Dynamic chromatin remodeling requires regulatory mechanisms for its adaptation to different nuclear function, which are likely to be mediated by signaling proteins. In this context, VRK1 is a nuclear Ser-Thr kinase that regulates pathways associated with transcription, replication, recombination, and DNA repair. Therefore, VRK1 is a potential regulatory, or coordinator, molecule in these processes. In this work we studied the effect that VRK1 depletion has on the basal nuclear and chromatin phosphoproteome, and their associated pathways. VRK1 depletion caused an alteration in the pattern of the nuclear phosphoproteome, which is mainly associated with nucleoproteins, ribonucleoproteins, RNA splicing and processing. Next, it was determined the changes in proteins associated with DNA damage that was induced by doxorubicin treatment. Doxorubicin alters the nuclear phosphoproteome affecting proteins implicated in DDR, including DSB repair proteins NBN and 53BP1, cellular response to stress and chromatin organization proteins. In VRK1-depleted cells, the effect of doxorubicin on protein phosphorylation was reverted to basal levels. The nuclear phosphoproteome patterns induced by doxorubicin are altered by VRK1 depletion, and is enriched in histone modification proteins and chromatin associated proteins. These results indicate that VRK1 plays a major role in processes requiring chromatin remodeling in its adaptation to different biological contexts.

Secretome processing for proteomics: A methods comparison.

The cancer cell secretome comprises a treasure-trove for biomarkers since it reflects cross-talk between tumor cells and their surrounding environment with high detectability in biofluids. In this study, we evaluated six secretome sample processing workflows coupled to single-shot mass spectrometry: (1) Protein concentration by ultrafiltration with a molecular weight cut-off (MWCO) filter and sample preparation through in-gel digestion (IGD); (2) Acetone protein precipitation coupled to IGD; (3) MWCO filter-based protein concentration followed by to in-solution digestion (ISD); (4) Acetone protein precipitation coupled to ISD; (5) Direct ISD; (6) Secretome lyophilization and ISD. To this end, we assessed workflow triplicates in terms of total number of protein identifications, unique identifications, reproducibility of protein identification and quantification and detectability of small proteins with important functions in cancer biology such as cytokines, chemokines, and growth factors. Our findings revealed that acetone protein precipitation coupled to ISD outperformed the other methods in terms of the number of identified proteins (2246) and method reproducibility (correlation coefficient between replicates (r = 0.94, CV = 19%). Overall, especially small proteins such as those from the classes mentioned above were better identified using ISD workflows. Concluding, herein we report that secretome protein precipitation coupled to ISD is the method of choice for high-throughput secretome proteomics via single shot nanoLC-MS/MS.

Candidate biomarkers for treatment benefit from sunitinib in patients with advanced renal cell carcinoma using mass spectrometry-based (phospho)proteomics.

The tyrosine kinase inhibitor sunitinib is an effective first-line treatment for patients with advanced renal cell carcinoma (RCC). Hypothesizing that a functional read-out by mass spectrometry-based (phospho, p-)proteomics will identify predictive biomarkers for treatment outcome of sunitinib, tumor tissues of 26 RCC patients were analyzed. Eight patients had primary resistant (RES) and 18 sensitive (SENS) RCC. A 78 phosphosite signature (p < 0.05, fold-change > 2) was identified; 22 p-sites were upregulated in RES (unique in RES: BCAR3, NOP58, EIF4A2, GDI1) and 56 in SENS (35 unique). EIF4A1/EIF4A2 were differentially expressed in RES at the (p-)proteome and, in an independent cohort, transcriptome level. Inferred kinase activity of MAPK3 (p = 0.026) and EGFR (p = 0.045) as determined by INKA was higher in SENS. Posttranslational modifications signature enrichment analysis showed that different p-site-centric signatures were enriched (p < 0.05), of which FGF1 and prolactin pathways in RES and, in SENS, vanadate and thrombin treatment pathways, were most significant. In conclusion, the RCC (phospho)proteome revealed differential p-sites and kinase activities associated with sunitinib resistance and sensitivity. Independent validation is warranted to develop an assay for upfront identification of patients who are intrinsically resistant to sunitinib.

Identification of protein biomarkers for prediction of response to platinum-based treatment regimens in patients with non-small cell lung cancer.

The majority of patients with resected stage II-IIIA non-small cell lung cancer (NSCLC) are treated with platinum-based adjuvant chemotherapy (ACT) in a one-size-fits-all approach. However, a significant number of patients do not derive clinical benefit, and no predictive patient selection biomarker is currently available. Using mass spectrometry-based proteomics, we have profiled tumour resection material of 2 independent, multi-centre cohorts of in total 67 patients with NSCLC who underwent ACT. Unsupervised cluster analysis of both cohorts revealed a poor response/survival sub-cluster composed of ~ 25% of the patients, that displayed a strong epithelial-mesenchymal transition signature and stromal phenotype. Beyond this stromal sub-population, we identified and validated platinum response prediction biomarker candidates involved in pathways relevant to the mechanism of action of platinum drugs, such as DNA damage repair, as well as less anticipated processes such as those related to the regulation of actin cytoskeleton. Integration with pre-clinical proteomics data supported a role for several of these candidate proteins in platinum response prediction. Validation of one of the candidates (HMGB1) in a third independent patient cohort using immunohistochemistry highlights the potential of translating these proteomics results to clinical practice.

Phosphoproteomics of patient-derived xenografts identifies targets and markers associated with sensitivity and resistance to EGFR blockade in colorectal cancer.

Epidermal growth factor receptor (EGFR) is a well-exploited therapeutic target in metastatic colorectal cancer (mCRC). Unfortunately, not all patients benefit from current EGFR inhibitors. Mass spectrometry-based proteomics and phosphoproteomics were performed on 30 genomically and pharmacologically characterized mCRC patient-derived xenografts (PDXs) to investigate the molecular basis of response to EGFR blockade and identify alternative drug targets to overcome resistance. Both the tyrosine and global phosphoproteome as well as the proteome harbored distinctive response signatures. We found that increased pathway activity related to mitogen-activated protein kinase (MAPK) inhibition and abundant tyrosine phosphorylation of cell junction proteins, such as CXADR and CLDN1/3, in sensitive tumors, whereas epithelial-mesenchymal transition and increased MAPK and AKT signaling were more prevalent in resistant tumors. Furthermore, the ranking of kinase activities in single samples confirmed the driver activity of ERBB2, EGFR, and MET in cetuximab-resistant tumors. This analysis also revealed high kinase activity of several members of the Src and ephrin kinase family in 2 CRC PDX models with genomically unexplained resistance. Inhibition of these hyperactive kinases, alone or in combination with cetuximab, resulted in growth inhibition of ex vivo PDX-derived organoids and in vivo PDXs. Together, these findings highlight the potential value of phosphoproteomics to improve our understanding of anti-EGFR treatment and response prediction in mCRC and bring to the forefront alternative drug targets in cetuximab-resistant tumors.

DPHL v.2: An updated and comprehensive DIA pan-human assay library for quantifying more than 14,000 proteins.

A comprehensive pan-human spectral library is critical for biomarker discovery using mass spectrometry (MS)-based proteomics. DPHL v.1, a previous pan-human library built from 1,096 data-dependent acquisition (DDA) MS data of 16 human tissue types, allows quantifying of 10,943 proteins. Here, we generated DPHL v.2 from 1,608 DDA-MS data. The data included 586 DDA-MS data acquired from 18 tissue types, while 1,022 files were derived from DPHL v.1. DPHL v.2 thus comprises data from 24 sample types, including several cancer types (lung, breast, kidney, and prostate cancer, among others). We generated four variants of DPHL v.2 to include semi-tryptic peptides and protein isoforms. DPHL v.2 was then applied to two colorectal cancer cohorts. The numbers of identified and significantly dysregulated proteins increased by at least 21.7% and 14.2%, respectively, compared with DPHL v.1. Our findings show that the increased human proteome coverage of DPHL v.2 provides larger pools of potential protein biomarkers.

Potential diagnostic value of CSF metabolism-related proteins across the Alzheimer's disease continuum.

BACKGROUND: Alzheimer's disease (AD) cerebrospinal fluid (CSF) core biomarkers (Aß42/40 ratio, p-tau, and t-tau) provide high diagnostic accuracy, even at the earliest stage of disease. However, these markers do not fully reflect the complex AD pathophysiology. Recent large scale CSF proteomic studies revealed several new AD candidate biomarkers related to metabolic pathways. In this study we measured the CSF levels of four metabolism-related proteins not directly linked to amyloid- and tau-pathways (i.e., pyruvate kinase, PKM; aldolase, ALDO; ubiquitin C-terminal hydrolase L1, UCHL1, and fatty acid-binding protein 3, FABP3) across the AD continuum. We aimed at validating the potential value of these proteins as new CSF biomarkers for AD and their possible involvement in AD pathogenesis, with specific interest on the preclinical phase of the disease. METHODS: CSF PKM and ALDO activities were measured with specific enzyme assays while UCHL1 and FABP3 levels were measured with immunoassays in a cohort of patients composed as follows: preclinical AD (pre-AD, n = 19, cognitively unimpaired), mild cognitive impairment due to AD (MCI-AD, n = 50), dementia due to AD (ADdem, n = 45), and patients with frontotemporal dementia (FTD, n = 37). Individuals with MCI not due to AD (MCI, n = 30) and subjective cognitive decline (SCD, n = 52) with negative CSF AD-profile, were enrolled as control groups. RESULTS: CSF UCHL1 and FABP3 levels, and PKM activity were significantly increased in AD patients, already at the pre-clinical stage. CSF PKM activity was also increased in FTD patients compared with control groups, being similar between AD and FTD patients. No difference was found in ALDO activity among the groups. UCHL1 showed good performance in discriminating early AD patients (pre-AD and MCI-AD) from controls (AUC ~ 0.83), as assessed by ROC analysis. Similar results were obtained for FABP3. Conversely, PKM provided the best performance when comparing FTD vs. MCI (AUC = 0.80). Combination of PKM, FABP3, and UCHL1 improved the diagnostic accuracy for the detection of patients within the AD continuum when compared with single biomarkers. CONCLUSIONS: Our study confirmed the potential role of UCHL1 and FABP3 as neurodegenerative biomarkers for AD. Furthermore, our results validated the increase of PKM activity in CSF of AD patients, already at the preclinical phase of the disease. Increased PKM activity was observed also in FTD patients, possibly underlining similar alterations in energy metabolism in AD and FTD.

Atrial epicardial adipose tissue abundantly secretes myeloperoxidase and activates atrial fibroblasts in patients with atrial fibrillation.

BACKGROUND: Epicardial adipose tissue (EAT) secretome induces fibrosis. Fibrosis, primarily extracellular matrix (ECM) produced by fibroblasts, creates a substrate for atrial fibrillation (AF). Whether the EAT secretome from patients with AF activates human atrial fibroblasts and through which components, remains unexplored. RESEARCH AIMS: (a) To investigate if the EAT secretome from patients with versus without AF increases ECM production in atrial fibroblasts. (b) To identify profibrotic proteins and processes in the EAT secretome and EAT from patients with, who will develop (future onset), and without AF. METHODS: Atrial EAT was obtainded during thoracoscopic ablation (AF, n = 20), or open-heart surgery (future onset and non-AF, n = 35). ECM gene expression of human atrial fibroblasts exposed to the EAT secretome and the proteomes of EAT secretome and EAT were assessed in patients with and without AF. Myeloperoxidase and neutrophil extracellular traps (NETs) were assessed immunohistochemically in patients with paroxysmal, persistent, future onset, and those who remain free of AF (non-AF). RESULTS: The expression of COL1A1 and FN1 in fibroblasts exposed to secretome from patients with AF was 3.7 and 4.7 times higher than in patients without AF (p < 0.05). Myeloperoxidase was the most increased protein in the EAT secretome and EAT from patients with versus without AF (FC 18.07 and 21.57, p < 0.005), as was the gene-set neutrophil degranulation. Immunohistochemically, myeloperoxidase was highest in persistent (FC 13.3, p < 0.0001) and increased in future onset AF (FC 2.4, p = 0.02) versus non-AF. Myeloperoxidase aggregated subepicardially and around fibrofatty infiltrates. NETs were increased in patients with persistent versus non-AF (p = 0.03). CONCLUSION: In AF, the EAT secretome induces ECM gene expression in atrial fibroblasts and contains abundant myeloperoxidase. EAT myeloperoxidase was increased prior to AF onset, and both myeloperoxidase and NETs were highest in persistent AF, highlighting the role of EAT neutrophils in the pathophysiology of AF.

Phosphoproteomics guides effective low-dose drug combinations against pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is a devastating disease with a limited set of known driver mutations but considerable cancer cell heterogeneity. Phosphoproteomics provides a readout of aberrant signaling and has the potential to identify new targets and guide treatment decisions. Using two-step sequential phosphopeptide enrichment, we generate a comprehensive phosphoproteome and proteome of nine PDAC cell lines, encompassing more than 20,000 phosphosites on 5,763 phospho-proteins, including 316 protein kinases. By using integrative inferred kinase activity (INKA) scoring, we identify multiple (parallel) activated kinases that are subsequently matched to kinase inhibitors. Compared with high-dose single-drug treatments, INKA-tailored low-dose 3-drug combinations against multiple targets demonstrate superior efficacy against PDAC cell lines, organoid cultures, and patient-derived xenografts. Overall, this approach is particularly more effective against the aggressive mesenchymal PDAC model compared with the epithelial model in both preclinical settings and may contribute to improved treatment outcomes in PDAC patients.

C1q is increased in cerebrospinal fluid-derived extracellular vesicles in Alzheimer's disease: A multi-cohort proteomics and immuno-assay validation study.

INTRODUCTION: Extracellular vesicles (EVs) may propagate and modulate Alzheimer's disease (AD) pathology. We aimed to comprehensively characterize the proteome of cerebrospinal fluid (CSF) EVs to identify proteins and pathways altered in AD. METHODS: CSF EVs were isolated by ultracentrifugation (Cohort 1) or Vn96 peptide (Cohort 2) from non-neurodegenerative controls (n = 15, 16) and AD patients (n = 22, 20, respectively). EVs were subjected to untargeted quantitative mass spectrometry-based proteomics. Results were validated by enzyme-linked immunosorbent assay (ELISA) in Cohorts 3 and 4, consisting of controls (n = 16, n = 43, (Cohort3, Cohort4)), and patients with AD (n = 24, n = 100). RESULTS: We found > 30 differentially expressed proteins in AD CSF EVs involved in immune-regulation. Increase of C1q levels in AD compared to non-demented controls was validated by ELISA (∼ 1.5 fold, p (Cohort 3) = 0.03, p (Cohort 4) = 0.005). DISCUSSION: EVs may be utilized as a potential biomarker and may play a so far unprecedented role in immune-regulation in AD.

Advancing wide implementation of precision oncology: A liquid nitrogen-free snap freezer preserves molecular profiles of biological samples.

PURPOSE: In precision oncology, tumor molecular profiles guide selection of therapy. Standardized snap freezing of tissue biospecimens is necessary to ensure reproducible, high-quality samples that preserve tumor biology for adequate molecular profiling. Quenching in liquid nitrogen (LN2 ) is the golden standard method, but LN2 has several limitations. We developed a LN2 -independent snap freezer with adjustable cold sink temperature. To benchmark this device against the golden standard, we compared molecular profiles of biospecimens. METHODS: Cancer cell lines and core needle normal tissue biopsies from five patients' liver resection specimens were used to compare mass spectrometry (MS)-based global phosphoproteomic and RNA sequencing profiles and RNA integrity obtained by both freezing methods. RESULTS: Unsupervised cluster analysis of phosphoproteomic and transcriptomic profiles of snap freezer versus LN2 -frozen K562 samples and liver biopsies showed no separation based on freezing method (with Pearson's r 0.96 (range 0.92-0.98) and >0.99 for K562 profiles, respectively), while samples with +2 h bench-time formed a separate cluster. RNA integrity was also similar for both snap freezing methods. Molecular profiles of liver biopsies were clearly identified per individual patient regardless of the applied freezing method. Two to 25 s freezing time variations did not induce profiling differences in HCT116 samples. CONCLUSION: The novel snap freezer preserves high-quality biospecimen and allows identification of individual patients' molecular profiles, while overcoming important limitations of the use of LN2 . This snap freezer may provide a useful tool in clinical cancer research and practice, enabling a wider implementation of (multi-)omics analyses for precision oncology.

A transformer architecture for retention time prediction in liquid chromatography mass spectrometry-based proteomics.

Accurate retention time (RT) prediction is important for spectral library-based analysis in data-independent acquisition mass spectrometry-based proteomics. The deep learning approach has demonstrated superior performance over traditional machine learning methods for this purpose. The transformer architecture is a recent development in deep learning that delivers state-of-the-art performance in many fields such as natural language processing, computer vision, and biology. We assess the performance of the transformer architecture for RT prediction using datasets from five deep learning models Prosit, DeepDIA, AutoRT, DeepPhospho, and AlphaPeptDeep. The experimental results on holdout datasets and independent datasets exhibit state-of-the-art performance of the transformer architecture. The software and evaluation datasets are publicly available for future development in the field.

Dysregulation of Mycobacterium marinum ESX-5 Secretion by Novel 1,2,4-oxadiazoles.

The ESX-5 secretion system is essential for the viability and virulence of slow-growing pathogenic mycobacterial species. In this study, we identified a 1,2,4-oxadiazole derivative as a putative effector of the ESX-5 secretion system. We confirmed that this 1,2,4-oxadiazole and several newly synthesized derivatives inhibited the ESX-5-dependent secretion of active lipase LipY by Mycobacterium marinum (M. marinum). Despite reduced lipase activity, we did not observe a defect in LipY secretion itself. Moreover, we found that several other ESX-5 substrates, especially the high molecular-weight PE_PGRS MMAR_5294, were even more abundantly secreted by M. marinum treated with several 1,2,4-oxadiazoles. Analysis of M. marinum grown in the presence of different oxadiazole derivatives revealed that the secretion of LipY and the induction of PE_PGRS secretion were, in fact, two independent phenotypes, as we were able to identify structural features in the compounds that specifically induced only one of these phenotypes. Whereas the three most potent 1,2,4-oxadiazoles displayed only a mild effect on the growth of M. marinum or M. tuberculosis in culture, these compounds significantly reduced bacterial burden in M. marinum-infected zebrafish models. In conclusion, we report a 1,2,4-oxadiazole scaffold that dysregulates ESX-5 protein secretion.

Integration of stool microbiota, proteome and amino acid profiles to discriminate patients with adenomas and colorectal cancer.

BACKGROUND: Screening for colorectal cancer (CRC) reduces its mortality but has limited sensitivity and specificity. Aims We aimed to explore potential biomarker panels for CRC and adenoma detection and to gain insight into the interaction between gut microbiota and human metabolism in the presence of these lesions. METHODS: This multicenter case-control cohort was performed between February 2016 and November 2019. Consecutive patients ≥18 years with a scheduled colonoscopy were asked to participate and divided into three age, gender, body-mass index and smoking status-matched subgroups: CRC (n = 12), adenomas (n = 21) and controls (n = 20). Participants collected fecal samples prior to bowel preparation on which proteome (LC-MS/MS), microbiota (16S rRNA profiling) and amino acid (HPLC) composition were assessed. Best predictive markers were combined to create diagnostic biomarker panels. Pearson correlation-based analysis on selected markers was performed to create networks of all platforms. RESULTS: Combining omics platforms provided new panels which outperformed hemoglobin in this cohort, currently used for screening (AUC 0.98, 0.95 and 0.87 for CRC vs controls, adenoma vs controls and CRC vs adenoma, respectively). Integration of data sets revealed markers associated with increased blood excretion, stress- and inflammatory responses and pointed toward downregulation of epithelial integrity. CONCLUSIONS: Integrating fecal microbiota, proteome and amino acids platforms provides for new biomarker panels that may improve noninvasive screening for adenomas and CRC, and may subsequently lead to lower incidence and mortality of colon cancer.

Inferring kinase activity from phosphoproteomic data: Tool comparison and recent applications.

Aberrant cellular signaling pathways are a hallmark of cancer and other diseases. One of the most important signaling mechanisms involves protein phosphorylation/dephosphorylation. Protein phosphorylation is catalyzed by protein kinases, and over 530 protein kinases have been identified in the human genome. Aberrant kinase activity is one of the drivers of tumorigenesis and cancer progression and results in altered phosphorylation abundance of downstream substrates. Upstream kinase activity can be inferred from the global collection of phosphorylated substrates. Mass spectrometry-based phosphoproteomic experiments nowadays routinely allow identification and quantitation of >10k phosphosites per biological sample. This substrate phosphorylation footprint can be used to infer upstream kinase activities using tools like Kinase Substrate Enrichment Analysis (KSEA), Posttranslational Modification Substrate Enrichment Analysis (PTM-SEA), and Integrative Inferred Kinase Activity Analysis (INKA). Since the topic of kinase activity inference is very active with many new approaches reported in the past 3 years, we would like to give an overview of the field. In this review, an inventory of kinase activity inference tools, their underlying algorithms, statistical frameworks, kinase-substrate databases, and user-friendliness is presented. The most widely-used tools are compared in-depth. Subsequently, recent applications of the tools are described focusing on clinical tissues and hematological samples. Two main application areas for kinase activity inference tools can be discerned. (1) Maximal biological insights can be obtained from large data sets with group comparisons using multiple complementary tools (e.g., PTM-SEA and KSEA or INKA). (2) In the oncology context where personalized treatment requires analysis of single samples, INKA for example, has emerged as tool that can prioritize actionable kinases for targeted inhibition.

Phosphoproteomic Analysis of FLCN Inactivation Highlights Differential Kinase Pathways and Regulatory TFEB Phosphoserines.

In Birt-Hogg-Dubé (BHD) syndrome, germline loss-of-function mutations in the Folliculin (FLCN) gene lead to an increased risk of renal cancer. To address how FLCN inactivation affects cellular kinase signaling pathways, we analyzed comprehensive phosphoproteomic profiles of FLCNPOS and FLCNNEG human renal tubular epithelial cells (RPTEC/TERT1). In total, 15,744 phosphorylated peptides were identified from 4329 phosphorylated proteins. INKA analysis revealed that FLCN loss alters the activity of numerous kinases, including tyrosine kinases EGFR, MET, and the Ephrin receptor subfamily (EPHA2 and EPHB1), as well their downstream targets MAPK1/3. Validation experiments in the BHD renal tumor cell line UOK257 confirmed that FLCN loss contributes to enhanced MAPK1/3 and downstream RPS6K1/3 signaling. The clinically available MAPK inhibitor Ulixertinib showed enhanced toxicity in FLCNNEG cells. Interestingly, FLCN inactivation induced the phosphorylation of PIK3CD (Tyr524) without altering the phosphorylation of canonical Akt1/Akt2/mTOR/EIF4EBP1 phosphosites. Also, we identified that FLCN inactivation resulted in dephosphorylation of TFEB Ser109, Ser114, and Ser122, which may be linked to increased oxidative stress levels in FLCNNEG cells. Together, our study highlights differential phosphorylation of specific kinases and substrates in FLCNNEG renal cells. This provides insight into BHD-associated renal tumorigenesis and may point to several novel candidates for targeted therapies.

Phosphoproteomic profiling of T cell acute lymphoblastic leukemia reveals targetable kinases and combination treatment strategies.

Protein kinase inhibitors are amongst the most successful cancer treatments, but targetable kinases activated by genomic abnormalities are rare in T cell acute lymphoblastic leukemia. Nevertheless, kinases can be activated in the absence of genetic defects. Thus, phosphoproteomics can provide information on pathway activation and signaling networks that offer opportunities for targeted therapy. Here, we describe a mass spectrometry-based global phosphoproteomic profiling of 11 T cell acute lymphoblastic leukemia cell lines to identify targetable kinases. We report a comprehensive dataset consisting of 21,000 phosphosites on 4,896 phosphoproteins, including 217 kinases. We identify active Src-family kinases signaling as well as active cyclin-dependent kinases. We validate putative targets for therapy ex vivo and identify potential combination treatments, such as the inhibition of the INSR/IGF-1R axis to increase the sensitivity to dasatinib treatment. Ex vivo validation of selected drug combinations using patient-derived xenografts provides a proof-of-concept for phosphoproteomics-guided design of personalized treatments.